Comment on "Diffusion Monte Carlo study of jellium surfaces: Electronic densities 

and pair correlation functions" 
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In a fixed-node diffusion Monte Carlo calculation of the total energy of jellium slabs, Acioli and 
Ceperley [Phys. Rev. B 54, 17199 (1996)] reported jellium surface energies that at low electron 
densities were significantly higher than those predicted in the local-density approximation (LDA) of 
density-functional theory. Assuming that the fixed-node error in the slab and the bulk calculations 
cancel out, we show that their data yield surface energies that are considerably closer to the LDA 
and in reasonable agreement with those obtained in the random-phase approximation. 

PACS numbers: 71.15.Mb, 71.45.Gm 



Acioli and Ceperlejii presented the results of fixed- 
node difussion Monte Carlo (DMC) calculations of the 
total energy of jellium slabs at five different electron 
densities. Assuming that the released-node correction is 
very small at the electron densities of interest, these au- 
thors extracted surface energies by substracting from the 
fixed-node slab energy the corresponding released-node 
bulk energies of Ceperley and Alder;^ as parametrized 
by Perdew and Zunger>2, They concluded that at low 
electron densities jellium surface energies are significantly 
higher than those predicted in the local-density approxi- 
mation (LDA) of density- functional theory (DFT). 

In this Comment, we show that combining fixed-node 
slab and release-node bulk energies results in substantial 
imprecision. Instead, we expect a large degree of can- 
cellation between the fixed-node error in the slab and 
the bulk calculations, and conclude that by substracting 
from the fixed-node slab energies the corresponding fixed- 
node bulk energies one obtains jellium surface energies 
that are considerably closer to the LDA^ and in reason- 
able agreement with those obtained in the random-phase 
approximation (RPA).^ 

Acioli and Ceperleji considered jellium slabs of den- 
sity riQ 3/47r(rsao)^ (ao is the Bohr radius), the 
thickness of the positive background in each slab being 
L = 7.21 (tsOo). Surface energies were then obtained 
from 
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where e^^""^ and e'''""^ represent slab and bulk energies per 
particle, respectively. Fixed-node energies per particle 
in bulk jellium were reported by Ceperleji^ to be higher 
than their release-node counterparts by 9 x 10~^ Ry for 
= 2 and 2 x 10^"* Ry for — 5. Hence, combining 
fixed-node slab and release-node bulk energies Eq. 
yields for =2 and 5 surface energies that are too large 
by ~ 150 and 5 erg/cm^, respectively. These fixed-node 
errors represent ^ 40% and ~ 20% of the LDA correla- 
tion energy for = 2 and 5, respectively. Furthermore, 



they are comparable to the deviation of the total surface 
energies reported in Ref. Q from those obtained in the 
LDA. 

In order to derive reliable surface energies from the 
fixed-node slab calculations of Ref. we require fixed- 
node energies of bulk jellium. Such calculations were 
reported by Ceperley^ and more recently by Kwon et ali 
for rs=l, 5, 10, and 20. The fixed-node Slater-Jastrow 
correlation energies of Kwon et ali can be parametrized 
in Perdew-Zunger'^ form, as foUows- 
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or in Perdew- Wang^ form, as follows 
^buik ^ _(o.12436 + 0.027404 r,) In (l + ^^^^ ) Ry, 



with 



(3) 



a = 7.5957 ry2+3.5876rs + 1.8207rf/2-f0.47746r2. (4) 

Both Eqs. l(2Jl and Q have been adjusted to fit the fixed- 
node correlation energies of Kwon et al7 for rs=l, 5, 
and 10. The fixed-node bulk energy e^""= = -0.01482 Ry 
obtained by Ceperley and Alder— for = 2.07 is under- 
estimated by Eq. ||2J) by 2 x 10~^ Ry and overestimated 
by Eq. @ by 2 x lO^^Ry. At r, = 20, Eqs. @ and 
P|l underestimate the fixed-node bulk energy of Kwon et 
ali by 2 X lO"'* Ry and 10"'' Ry, respectively. 

Adding to the correlation energy of Eqs. |2Jl and 
the well-known kinetic and exchange energies of a uni- 
form electron gasiS and combining these energies with the 
fixed-node slab energies of Acioli and Ceperley,i Eq. 
yields the surface energies (ai and (72) shown in Table I. 
Also shown in this Table are the surface energies re- 
ported in Ref. I* {(Jac), together with the LDA and RPA 
surface energies that we have obtained for jellium slabs 
with L = 7.21 (ts ao). Small differences of no more than 
10 erg/cm^ between these LDA and RPA surface energies 



TABLE I: Surface energies of jellium slabs with L = 
7.21 (r^ ao), as obtained from Eq. Q by combining the Acioli- 
Ceperley fixed-node slab energies with the fixed-node bulk 
energies of either Eq. Q {oi) or Eq. Q (0-2), and by combin- 
ing the Acioli-Ceperley fixed-node slab energies with release- 
node bulk energies {(tac)- At = 2.07, combining the 
Acioli-Ceperley fixed-node slab energy with the correspond- 
ing fixed-node bulk energy = -0.01482 Ry of Ceperley 
and Alder^^ yields a = —558. ctlda and ctrpa represent LDA 
and RPA surface energies, ctq = ""s + o"es + represents the 
combined kinetic [a a), electrostatic (tJes), and exchange (era;) 
contributions to the total surface energy a. The correlation 
surface energy is simply CTc ~ a — oq. Units are erg/cm^ 
(1 erg/cm^ = 6.2415 x 10~^ eV/A^). 
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and those reported before for the semi-infinite jellium^ 
are due to the finite size of our jeUium slabs. 

An inspection of Table I shows that using in Eq. 
the fixed-node bulk energies of either Eq. |(2Jl or Eq. Q 
brings the DMC surface energies closer to the LDA and 
to reasonable agreement with the RPA. This is consistent 
with recent work, where it was shown that upon surface 
formation there is a persistent cancellation of short-range 
correlation effects beyond the RPAiii Other approaches 
have also led to the conclusion that the actual jellium 
surface energies should be only slightly higher than those 
obtained in the LDAia*iii4 

Li et al}^ calculated the fixed- node DMC surface en- 
ergy of a jellium slab with — 2.07, the thickness of 
the positive background being L — 8.52 (rsOo). These 
authorsr^ unlike Acioli and Ceperley, "'^ extracted from 
their fixed-node slab energy the corresponding fixed- 
node bulk energy (e''""= = -0.01482 Ry)V& and found 
a = — 465 erg/cm^jii they also performed LDA calcula- 
tions and found an LDA surface energy of —567 erg/cm^. 
We have calculated the LDA and RPA surface energies 
of a jelHum slab with Ts = 2.07 and L = 8.52 (r^ao), 
and have found ctrpa = —485 erg/cm^ and ctlda — 
—589 erg/cm^ jASt which are both only ^ 20 erg/cm^ 
smaller than the corresponding DMC surface energy re- 
ported in Ref. 0. 

Fig. 1 exhibits the surface energies of Table I, as a 
function of Vs. At the highest densities the agreement 
between DMC and RPA surface energies is good. At the 
smallest densities {vg > 3.25) the corrected DMC sur- 
face energies are still larger than their RPA counterparts. 
Nevertheless, we note that DMC surface energies are very 
sensitive to small uncertainties in the parametrization of 
the bulk energy. Therefore, if one is to quantitatively 
account for jellium surface energies, both DMC slab and 
bulk energies entering Eq. |^ should be computed on the 
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FIG. 1: Surface energies of jellium slabs with L — 7.21 (r^ ao), 
as obtained from Eq. Q. Stars and open circles repre- 
sent DMC surface energies obtained by combining the Acioli- 
Ceperley fixed-node slab energies with the fixed-node bulk 
energies of Eqs. Q and 0, respectively. The solid circle 
represents the DMC surface energy obtained at rs = 2.07 by 
combining the Acioli-Ceperley fixed-node slab energy with the 
corresponding fixed-node bulk energy e'"'"° = —0.01482 Ry of 
Ceperley and Alder^ Crosses represent the surface energies 
obtained in Ref. Q by combining fixed-node slab and release- 
node bulk energies. Dashed and solid lines represent LDA 
and RPA calculations, respectively. 
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FIG. 2: Correlation contribution to the surface energy. Stars 
and open circles represent DMC correlation surface energies 
obtained by combining the Acioli-Ceperley fixed-node slab en- 
ergies with the fixed-node bulk energies of Eqs. (|5J and lO, 
respectively. The solid circle represents the DMC correla- 
tion surface energy obtained at Vs = 2.07 by combining the 
Acioli-Ceperley fixed-node slab energy with the correspond- 
ing fixed-node bulk energy = -0.01482 Ry of Ceperley 
and Alder. Dashed and solid lines represent LDA and RPA 
calculations, respectively. 
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same footing. 

Small differences between the DMC and RPA calcula- 
tions come from the correlation contribution to the sur- 
face energy. We have carried out calculations of exact 
kinetic (cs), electrostatic (ues), and exchange (a^) sur- 
face energies for jellium slabs with L = 7.21 (r^ ao), and 
we have defined the correlation contribution to the DMC 
surface energy as 

(Tc ^ <T - {as + (Tes + <Jx)- (5) 

The results we have obtained from our analysis of the 
Acioli-Ceperley DMC slab energies are shown in Fig. 2 
(see also Table I), together with the LDA and RPA corre- 
lation energies that we have obtained for the same jellium 
slabs. This figure clearly shows that combining fixed- 
node slab energies with their fixed-node bulk counter- 
parts yields DMC correlation surface energies in reason- 
able agreement with the results obtained in the RPA. 
While the LDA underestimates considerably the correla- 
tion surface energy for all electron densities, it overesti- 
mates the exact exchange surface energy. Large and op- 
posite separate corrections to the LDA for exchange and 
correlation largely compensate, as discussed before^*^ 



and yield LDA surface energies that are close to their 
non-local counterparts. 

In conclusion, assuming that the fixed-node error in the 
slab and the bulk calculations cancel out, the DMC data 
reported in Ref. yields surface energies that are con- 
siderably closer to the LDA and in reasonable agreement 
with those obtained in the RPA. Nevertheless, at the 
smallest densities the corrected DMC surface energies are 
still larger than their RPA counterparts; at these densi- 
ties, they are also larger than the jellium surface energies 
extracted from DMC calculations for jellium spheres^ 
which are found to be close to the LDAii4 Surface ener- 
gies are extremely sensitive to little uncertainties in both 
slab and bulk energies; hence, in order to quantitatively 
account for the impact of nonlocal xc effects on the sur- 
face energy one needs to pursue DMC evaluations of both 
slab and bulk energies on the same footing. 
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